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In this paper, we calculate the spectrum of scalar field fluctuations in a bouncing, asymptotically 
flat Universe, and investigate the dependence of the result on changes in the physics on length scales 
shorter than the Planck length which are introduced via modifications of the dispersion relation. In 
this model, there are no ambiguities concerning the choice of the initial vacuum state. We study an 
example in which the final spectrum of fluctuations depends sensitively on the modifications of the 
dispersion relation without needing to invoke complex frequencies. Changes in the amplitude and 
in the spectral index are possible, in addition to modulations of the spectrum. This strengthens 
the conclusions of previous work in which the spectrum of cosmological perturbations in expanding 
inflationary cosmologies was studied, and it was found that, for dispersion relations for which the 
evolution is not adiabatic, the spectrum changes from the standard prediction of scale-invariance. 

PACS numbers: 98.80.Cq, 98.70.Vc 


I. INTRODUCTION 


The dependence of the spectrum of cosmological fluc¬ 
tuations in an inflationary Universe on hidden assump¬ 
tions about the physics on length scales much smaller 
than the Planck length has been recently studied in Refs. 
10, 0, H]. A priori this dependence comes about since in 
typical scalar-field-driven inflationary models the dura¬ 
tion of the period of expansion is so long that the phys¬ 
ical wavelengths of comoving modes which correspond 
to the present large-scale structure of the Universe are 
much smaller than the Planck length at the beginning 
of inflation. In weakly coupled scalar field models of in¬ 
flation, the spectrum of fluctuations is calculated by as¬ 
suming that fluctuation modes start out in the vacuum 
(i.e. minimum energy density) state at the beginning of 
inflation and subsequently evolve as determined by the 
equations of motion. It is hence not unreasonable to ex¬ 
pect that modihcations of the physics on trans-Planckian 
scales could affect the final spectrum of fluctuations. 

While the correct theory of trans-Planckian physics is 
not known, possible effects of the new physics can be 
modeled by changes in the dispersion relation of the helds 
corresponding to linear cosmolomcal fluctuations. Such 
an approach was used earlier 0, 0] to study the possible 
dependence of the spectrum of black hole radiation on 
trans-Planckian physics. Indeed, for a class of dispersion 
relations which deviate so strongly from the usual one 
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such that the evolution of the mode functions is highly 
non-adiabatic and the effective frequency becomes imag¬ 
inary, it was found that the spectrum of fluctuations is 
modified in a signihcant way compared to what is ob¬ 
tained for a linear dispersion relation 0 0 . As shown in 
0, 1^ , the final spectrum of fluctuations is unchanged 

if the change in the dispersion relation is not too drastic, 
the wave function evolves adiabatically, and hence the 
WKB approximation is valid throughout the evolution. 

Two major deficiencies in the previous work concerned 
the choice of initial conditions and the fact that com¬ 
plex dispersion relations were required in order to obtain 
signihcant deviations from the standard results. In the 
previous papers, the assumption was made that the state 
starts out as the one minimizing the energy density. How¬ 
ever, if the dispersion relation differs greatly from the 
linear one, in particular if it becomes complex, the phys¬ 
ical motivation for choosing this state becomes unclear. 
This problem does not arise in a bouncing Universe. In 
this case, modes corresponding to present cosmological 
scales could well have had a physical wavelength smaller 
than the Planck length at the bounce point, but at early 
(pre-bounce) times, the wavelength was larger than the 
Planck length. In this context, initial conditions can be 
set up at some very early pre-bounce time during which 
the mode obeys a linear dispersion relation and the physi¬ 
cal meaning of the minimum energy density state is clear. 
It is thus interesting to study the spectrum of fluctua¬ 
tions in a bouncing Universe which is asymptotically flat 
such that the choice of initial vacuum becomes unam¬ 
biguous. It is also interesting to exhibit a case where 
complex dispersion relations are not needed. Although 
complex frequencies are standard in classical physics and 
in quantum mechanics, they seem to be more problematic 
in the context of quantum field theory. 

In this article, we find an example of a bouncing Uni¬ 
verse in which the spectrum of fluctuations depends sen- 
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sitively on changes in the physics at length scales smaller 
than the Planck length, without requiring complex fre¬ 
quencies at high wavenumbers, thus improving on the 
second major deficiency of previous work. Another rea¬ 
son for investigating the dependence of the spectrum of 
fluctuations on modifications of trans-Planckian physics 
in a bouncing Universe is that one could expect (how¬ 
ever, as we show, this is not correct) that the devia¬ 
tions in the spectrum incurred during the pre-bounce 
period when the wavelength is blue-shifting in the trans- 
Planckian regime will cancel with the deviation in the 
post-bounce trans-Planckian regime when the mode is 
being red-shifted. Such a result would then lead to 
the expectation that the results on the dependence of 
the spectrum of fluctuations on changes in the trans- 
Planckian physics might be different in the case of black 
holes and (non-bounce) inflationary cosmology. A third 
reason for analyzing a bouncing Universe is that there 
has recently been a lot of interest in the field of string 
cosmology in bouncing Universe models. In the Einstein 
frame, the pre-big-bang scenario [|^ starts in a collaps¬ 
ing dilaton-dominated phase, and the same is true in the 
Ekpyrotic scenario |ll[| . Our arguments show that trans- 
Planckian effects could change the predictions of stan¬ 
dard cosmological perturbation theory in these examples 
1 


A model for the evolution of the scale factor in a bounc¬ 
ing inflationary Universe is a{t) = aocosh(ijf), where oq 
denotes the value of a{t) at the center of the bounce, 
and H is the Hubble expansion (contraction) rate dur¬ 
ing the period of exponential expansion (contraction). 
The idea of the calculation is to assume an initial spec¬ 
trum of fluctuations at some early time —t{ when the 
wavelength of all modes of interest is much larger than 
the Planck length, to propagate the perturbation modes 
through the bounce and calculate the final spectrum at 
time t{. In this model, one of the difficulties encountered 
in the previous work |Q, namely the problem of spec¬ 
ifying the initial state in the trans-Planckian regime, is 
overcome. However, it is necessary now to set up initial 
conditions on wavelengths which typically also are much 
larger than the Hubble radius. The latter problem can 
be overcome by considering a bouncing Universe which is 
radiation-dominated in the asymptotic past and future, 
e.g. 0,0 a{r]) = o? + where a and /3 are 

constants. Another toy model in which it is possible to 
set up well motivated initial conditions on scales much 
smaller than the Hubble radius and at the same time 
much larger than the Planck scale is a model in which 
the Universe is asymptotically flat both at large positive 


and negative times: 


a(ry) 


4 - 


^0 — 4 

1 + {v/m? 


( 1 ) 


where Iq and 4 (with £b (.q) denote the asymptotic 
size of the Universe and the size at the bounce, respec¬ 
tively, and 770 determines the time scale of the bounce, 
i.e. the time over which the scale factor changes. In 
the above, we are using conformal time rj, and employing 
the convention that the scale factor carries dimension of 
length while the comoving coordinates (including confor¬ 
mal time) are dimensionless. 

Finally, a comment is in order on the issue of backre- 
action. It has been shown in Refs. i 0 that the en¬ 
ergy density of the particles created by trans-Planckian 
physics can contribute in a significant way to the back¬ 
ground energy density. We do not address this interesting 
question since this clearly comes as a second step once 
a case where the spectrum is modified has been explic¬ 
itly exhibited. In this paper we concentrate on the first 
question since it would be useless to treat the backre- 
action question if no well-motivated example where the 
spectrum is modified can be found. However, it should 
be clear that, once this is done (as shown in the present 
article), this problem becomes the central question in the 
study of trans-Planckian physics in cosmology. 

This article is organized as follows. In section H, we 
review the arguments which show that no change in the 
final spectrum of fluctuations of a scalar matter field 
in a bouncing Universe is expected independent of the 
dispersion relation, provided that the WKB approxima¬ 
tion remains justified. We give a qualitative reason to 
expect changes in the amplitude of the final spectrum 
for dispersion relations for which the adiabatic approxi¬ 
mation breaks down. As mentioned above, our analysis 
assumes that perturbation theoryremains justified and 
back-reaction can be neglected [0. In Section HI we 
study a concrete model, namely the asymptotically flat 
bouncing Universe given by (^ with a dispersion relation 
modified in the trans-Planckian regime according to the 
prescription of Corley and Jacobson |^. We show that, 
as expected, both the overall amplitude and the shape of 
the spectrum differ from what is obtained for the linear 
dispersion relation, without ever requiring the effective 
frequency to become imaginary. The change in the spec¬ 
trum is produced by an interesting combination of effects 
due to the modified dispersion relation and the driving 
term for Parker particle production [p0| . 


II. METHOD AND QUALITATIVE 
CONSIDERATIONS 


^ Note that there is disagreement on the result of the linear theory The equation of motion for a minimally COUpled free 

of fluctuations in the Ekpyrotic scenario ^ scalar matter field ^{rj ,) in a Universe with scale factor 
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a(r]) in momentum space takes on the simple form 

= 0 , ( 2 ) 

making use of the re-scaling 

J (3) 

where n denotes the comoving wavenumber linked to the 
physical wavenumber by the relation k = n/a{ri). In 
the above, a prime denotes the derivative with respect 
to conformal time. The previous equations are valid for 
a spacetime with flat spatial sections. The advantage of 
considering a test scalar field and/or gravitational waves 
is that we do not need to address the origin of the dynam¬ 
ics of the scale factor. This is in contrast to the case of 
density perturbations. However, the Friedmann equation 
reads H?jc? -I- Kja} = where H. = a!ja and 

Pi > 0 is the energy density. This means that a bounce 
is consistent with the Einstein equations only il K = 1. 
For iF = 0, it would imply p = 0 for each component. 
Therefore, it seems that it is inconsistent to use Eq. (H) 
if the scale factor behaves such that there is a bounce. 
There are two ways out of the previous argument. The 
first one is the following. UK = 1, then Eq. (^) takes 
the form p" -|- {v? — K — a"/a)fi = 0. In addition the 
eigenfunctions of the three-dimensional Laplacian oper¬ 
ator are no longer planes waves. However, if we restrict 
ourselves to modes with wavelengths much smaller than 

I 


the curvature radius, then they do not feel the curva¬ 
ture of the spacelike section and we can safely work with 
Eq. (^). It could be checked that this is indeed the case 
for the wavenumbers considered in the example presented 
in the following section. Another possibility is to use the 
“nonsingular Universe construction” [^. In this case, it 
is definitely possible to get a spatially flat bounce. What 
happens is that the higher derivative terms become im¬ 
portant at the bounce and enable the evasion of the pre¬ 
vious argument by supplying other terms in the analog of 
the first Friedmann equation. In conclusion, although in 
this article we consider a scale factor with a bounce, we 
can nevertheless work with Eq. (^, i.e. we let the non¬ 
singular Universe construction enter only via the scale 
factor. 

The method introduced in i§ to study the depen¬ 
dence of the spectrum of fluctuations on trans-Planckian 
physics is to replace the linear dispersion relation t'-’phys = 
A: by a non standard dispersion relation (fc) 

where the function is a priori arbitrary. In the 

context of cosmology, it has been shown in Ref. |Q, ||] 
that this amounts to replacing appearing in (||) with 
nlff{n,ri) defined by 




a{r]) 


( 4 ) 


For a fixed comoving mode, this implies that the disper¬ 
sion relation becomes time-dependent. Various forms for 
rieff have been used, e.g. 




n^tanh^/P 



p-l 


the first being a generalization of the relation studied by 
Unruh [|j whereas the second one is a generalization of 
the relations used by Corley and Jacobson Q. In the 
above, \{Tf) is the physical wavelength, is the cutoff 
length (which can be taken to be the Planck length), p 
and m are integers, and bm is a real number which can be 
either positive or negative. Thus, the equation of motion 
to be analyzed is 


nlaiji.r]) - 


fj, = 0. 


( 6 ) 


For fixed comoving wavenumber n, the evolution of /i 
depends crucially on which region rj lies in. The hrst 
region is the trans-Planckian region with 


X{r]) < 4 (Region I). (7) 


Let us first consider the standard inflationary scale factor 
and a monotonic dispersion relation [like the first one in 


n-lsin,r)) = 




Xv). 


( 5 ) 


Eq. (||) or the second one with > 0]. This means 
that, initially, the term ngjj(n, 77 ) dominates in Eq. ®)- 
The initial conditions are fixed in this region and since 
the WKB approximation is applicable, we can choose the 
initial state as the “minimizing energy state” [^ . Then, 
the (positive frequency) solution is given by 

, ^ 1 

n,v) 

where r]i is some initial time. 

The second region corresponds to 

X ^ and a"/a (Region 2). (9) 

In this region the mode has reached the linear part of the 
dispersion relation. The general solution in Region 2 is 
the plane wave 

pL = Bi exp{—inr]) + B 2 exp(zn77) (10) 


H-vil — 


\J 2 neff ( 


exp 


—i 


Tiesdr 


( 8 ) 
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with constant coefficients Bi and -62 • For the stan¬ 
dard dispersion relation, the initial conditions are fixed 
in this region. The usual choice of the vacuum state is 
Bi = 1/-^/^, B 2 = 0. In general Bi and B 2 are de¬ 
termined by the matching conditions between Regions 1 
and 2. However, if the dynamics is adiabatic throughout 
(in particular if the a"/a term is negligible), the WKB 
approximation holds and the solution is always given by 
(^. Therefore, if we start with B 2 = 0 and uses this so¬ 
lution, one finds that B 2 remains zero at all times. Deep 
in the region where riefi ~ n the solution becomes 




—f= ex.p{-i(j) 
y2n 


inrf)^ 


( 11 ) 


i.e. the standard vacuum solution times a phase which 
will disappear when we calculate the modulus. The phase 
(j) is given by 0 = Ueffdr, where rji is the time at which 
rieff ~ n. 

The situation changes dramatically if we consider non¬ 
monotonic dispersion relations. This is the case if 6 ^ < 0 
in Eq. (|^) or for the dispersion relation introduced in Ref. 
[p^ . Two new features can occur. Either the dispersion 
relation can become complex (this is in general the case 
if bm < 0 ) or the term a”/a dominates (this is the case 
for the dispersion relation introduced in Ref. or for 
the Cor ley/Jacobson dispersion relation with bm < 0 for 
a certain range of comoving wavenumbers). In this case, 
the WKB approximation is violated in Region 1 and we 
expect changes in the final spectrum. Unfortunately, in 
the context of the standard increasing inflationary scale 
factor, one also looses the ability to fix natural initial con¬ 
ditions. However, this is no longer true if the spacetime 
is asymptotically flat because then, at infinity, the term 
a”/a goes to zero, see also Ref. Q. Therefore, in this 
case, we can choose well-motivated initial conditions. A 
good example of such a situation is provided by a bounc¬ 
ing Universe. From the above considerations, we expect 
that in this case the final spectrum is modified and the 
initial conditions can be fixed naturally. As an additional 
benefit, complex frequencies can be avoided. In the next 
section, we consider a toy model where these qualitative 
arguments can be implemented concretely, at the level of 
equations. 


III. A SPECIFIC EXAMPLE 

We will now illustrate the qualitative arguments of the 
previous section with a concrete quantitative example. 
We take the asymptotically flat bouncing Universe given 
by d^). We consider the type of modified dispersion re¬ 
lation for which in the case of an expanding inflation¬ 
ary Universe the deviations in the spectrum were found 
[|^, ^, namely the (generalization of the) Cor ley/Jacobson 
dispersion relation (||) with b^ negative. In this case, 
Wphys(fc) is given by 

= > ( 12 ) 

where is the cutoff physical wavenumber. The disper¬ 
sion relation is represented in Fig. |^. 



FIG. 1: Corley/Jacobson dispersion relation for m = 1 and 
\bm\ = 1. In general Wpj^y^(fc) vanishes at fc = ■ 

For k > the physical frequency becomes imag¬ 

inary. 

The behavior of the solutions of the equation of motion 
/i" + — a" I a) = 0 is determined by the competition 

between the two terms n^g and a”/a given by 


J 


nlsin, T]) 


uV" 

-1 

■Si 

0 

1 

cr 

—2m 

a" _ 2(4 - 4) 1 - 3(77/770)^ 


r i+{v/vor\ 

5 

a vUh [1 + {ri/mY? 


(13) 


where = 1/k^ is the cutoff length. These two terms are represented in Fig. ||. When n^g > a"/a, which is always 
the case when r]/r]o —> ± 00 , the WKB approximation is valid and the fundamental solutions can be written as 


fj. ~ 




exp 


±i 


neff(n,T)dT 

I— 


^'n/rio^±oc 


^/2n 


^±inri 


(14) 


In this case, the physical wavenumbers of the modes are such that they correspond to the linear part of the disper- 
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sion relation, see Fig. |^. In this situation, natural initial 
conditions can be chosen. Note that these initial condi¬ 
tions are, in a sense, even “more standard” than in the 
previous studies on trans-Planckian physics since usually 
the initial conditions are set in a region where the disper¬ 
sion relation is not linear (but, as explained above, since 
the WKB approximation holds, meaningful initial condi¬ 
tions can nevertheless be considered). On the contrary, 
when rigfj < a!' ja the WKB approximation is violated 
and two independent solutions are given by 

dT 

^-0(77)/ ( 15 ) 

J a^(T) 


In this case, we are close to the cutoff scale and (k) 
being close to zero, see Fig. the a"/a term dominates. 
This situation is very similar to the one discussed in Ref. 

where this part of the dispersion relation has been 
named “the tail”. In fact, it is possible to quantify ex¬ 
actly the accuracy of the WKB approximation. Given an 
equation of the form /r" -|- = 0, the WKB approx¬ 
imation is valid if the quantity 1, where the 

quantity Q is defined by the following expression 


_ 3(w')" 

4^2 2lo ■ 


(16) 


This standard criterion can be obtained in the following 
manner. The WKB solution, /Uwkb, satisfies the equation 
Mwkb + Mwkb(i^^ — Q) = 0 exactly. Therefore, one has 
/r ~ /Twkb if IQ/w^l <C 1. In the present context, is of 
course equal to — a"/a. We have plotted the quantity 
\Q/ijj"^\ in Fig. When |77/?7o| > 1, < 1 and the 

WKB approximation is satisfied, as expected from the 
previous discussion. The quantity \Q/uP'\ increases when 
one approaches the bounce. At rj = rjj, where = 

0, there is a divergence. This point is a turning point. 
Then for 177 / 770 ! 1) as can be seen in Fig. ||, jQ/w^j 

remains large and the WKB approximation breaks down 
in agreement with the considerations above. 

In a bouncing Universe, the modes start in the linear 
part of the dispersion relation, pass through the non¬ 
linear part and then come back in the standard region. 
However, things are not so simple because there exists a 
range of comoving wavenumbers such that the dispersion 
relation becomes complex. Although this case is a priori 
interesting, it clearly requires more speculative consider¬ 
ations, in particular the quantization in the presence of 
imaginary frequency modes. Since our goal in this pa¬ 
per is to exhibit a case where everything can be done in 
a standard manner, we will restrict ourselves to modes 
which never enter the region where the dispersion rela¬ 
tion becomes complex. In the following, we determine the 
range of comoving wavenumbers we are interested in and 
for which we are going to calculate the power spectrum. 

The maximum of the absolute value of a"/a is located 
at 77 = 0 and is given by ip jj — 1). Therefore, if 





FIG. 2: The solid curve is a!' j a, the dashed one Weff. The 
values used are Iq = 5000, = 10, = 1, jbmj = 1, tti = 1 

and n — 60. 



FIG. 3: Absolute value of Q versus conformal time. The 
values used are = 5000, Ih = 10, = 1, jbmj = 1, m = 1 

and n = 55. 


we restrict ourselves to modes such that 


n :$> Uh = 


Vo \4 


1 


(17) 


then one is sure that, with an unmodified dispersion re¬ 
lation, the term a"/a can always be neglected and that 
the initial spectrum is never changed. Clearly, this is not 
a physical restriction but it renders the comparison with 
the case with a modified dispersion relation easier. The 
minimum of the modified part of the dispersion relation is 
given by (from now on, we consider the case \bm\ = 1 and 
777 = 1 since it is more convenient and does not restrict 
the physical content of the problem in any way) 


nls{n,V = 0) 



Uinf = 27r^. (18) 


Therefore, to maintain a real dispersion relation, we 
should only consider modes such that n < riinf. Of 
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course, for consistency, the parameters of the model must 
be chosen such that Ub < Uinf. This is the case for the 
most natural choice, i.e., ^ ^ -^c- the following, 

as already mentioned above, the time jyj (n) such that 


that the range reduces to 

nl<n^<nt,,. ( 22 ) 


nls{n,Vi) = ^ ( 19 ) 

will play a crucial role. For convenience we will only 
consider values of n such that this time is determined only 
by the central peak of a”/a and not by the two wings. In 
practice, this amounts to taking oo) > rig/A since 

n^/4 is the maximum height of the wings. Then it is easy 
to show that G [n?_, where 


Having determined the relevant wavenumbers, we can 
now choose the initial condition and solve the equation 
of motion. In the first region where rj < —rj^(n), we only 
consider positive frequency modes and we have 


A^i(^) = 




exp 


rv 

i / neff(n,r)dr 


(23) 



where Ugup = 27r£o/^c 1 ® largest value of n such that 
ngg(n, oo) remains real. Therefore, the range we are in¬ 
terested in is given by 

max(nb,n.?_) < < min(n?„f, n^). ( 21 ) 

In practice, we have n‘/_ ~ (I/4)n^ and n\ ~ ^T'sup- This 
means that max(n^,n?_) = rig, min(n?jf,n^) = and 


where rji is an arbitrary initial time. In the second region, 
where —r]j{n) < rj < rij(ji), the solution is given by 

miv) = Bia{r]) + B2a{r]) [ (24) 

Jo a (t) 

The lower bound of the integral is a priori arbitrary. 
However, it is very convenient to take it equal to zero be¬ 
cause in this case the second branch becomes odd whereas 
the first one (i.e. the scale factor) is even. Then, it is 
easy to show that 


r]o( 


-b 


1 — 2p + p^ 


n dr 

/o (B{t) (IXp'^Po ' 2p'^[l+p{r)/'qoY]r]^ 


T] P^ +2p- 3 f T] 

+ 2 p ./2 


where p = Iq/Io- Finally, the solution in the third region where rj > rij{n) can be written as 


foii(’7) = 


Cl 


'\/2?T.efF 


exp 


-i / neff(n,T)dT 

^r]i 


c. 




exp 


+i / neff(«,'r)dT 
J r/i 


(25) 


( 26 ) 


The goal is now to calculate the coefficients Ci and C 2 - Using the continuity of the mode function /r and of its 
derivative, we find 


^ md-rri / . ^ ^ + “S'](-??j)[/' + afj(vj) + [/' -b af](-rjj)[g' + agj(pj) L (27) 

2Wl-gj(n)j y'nefi(n,-?7j)neff(n, ?7j) L J 

^ md-rri / . ^ ^ T1 + “S'](-??j)[/' + a/](»?j) + [/' + afK-Vj)lff' + ag](rjj) I, (28) 

2M/[-77j(njJ y'nefi(n,-?7j)neff(n, ?7j) ( J 

where we have used the short-hand notation / = a(g) and g = 0 ( 77 ) Jg dT/a^(T). We have also utilized the following 
definitions: W = gf — g'f is the Wronskian, ^{rj) = nefi(n,T)dr, and the quantity a is 

a = + iries. (29) 

Using the fact that ries is even, it is easy to see that a{n,—g) = —a{n,r]). Using the same property, one could also 
simplify the factor —r]j)neg{n, pj) = nes{n,gj) in the above equation. The final result is given by Eqs. (|^), 

( p^ where all the functions are explicitly known except the function rjj = gj(ri). The dependence on the dispersion 
relation of the final result is completely encoded in this function. The time pj (n) is the solution of an algebraic 
equation that is not possible to solve explicitly. However, we can find an approximation if the function neff(n,, ry) is 
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Taylor expanded around rj = 0. For the function a"/a , we use a quadratic least square approximation which gives a 
better result. We obtain 




n- 


riinf J 


(l-«b’7 ), — - 


l- 6 (^ 

\Vo 


2 n 


(30) 


The two functions and their approximations are represented in Fig. We see now that ? 7 j(n) is the solution of a 
quadratic algebraic equation. Straightforward calculations give 


^•b L V J 


f jlV+ f + F 


n ) 


?7oJ 


(31) 


The function rjj (n) can also be found exactly by numerical 
calculations. We have done this computation by means 



V/Vo 

FIG. 4: The solid curve represents a" ja, the dot-dashed curve 
its quadratic approximation. The dashed curve represents 
and its quadratic approximation is the fourth curve. Val¬ 
ues chosen are lo = 5000, Ih = 10, = 1, |6m| = 1, m = 1 

and n — 60. 


of a Fortran code. The comparison of the exact result 
with its approximation given in Eq. (|^) is represented 
in Fig. ||. The function neff{n,r] = 0) is not monotonic 
and has a maximum around n ~ 45 for the values of the 
parameters considered in Fig. For this value nes{n,ri) 
is almost never smaller than a" j a. It is clear that around 
this value the approximation will be very good whereas 
for other values of n, especially for n ~ nb, the approxi¬ 
mation will be less good. This is due to the fact that when 
n approaches rib, the asymptotic value of nes{ri,r]) de¬ 
creases and the curve opens out at the intersection with 
a"/a. As a consequence, the quadratic approximation 
employed above breaks down. This is confirmed by the 
plots in Fig. |^. The approximation can be less good for 
other values of the parameters. 


From the above analysis, it is clear that e = rij{n )/tjq 
is a small number. We can therefore expand everything 
in terms of the parameter e. This will allow us to get an 
analytical estimate of the spectrum. After tedious but 
straightforward computations, one obtains the following 
expressions for the coefficients Ci{n) and (72(n) 


J 


Ci{n) 


C2{n) 


2n^ 1 




2n 


H? \ 


'-inf/ 


2 n^ 1 - 


1/2 



-1/2 

+ «b^ 



1/2 



«?nfy 

) 

-h nj; 1^1 — 

<1. 


1 - 


'-inf 


- 1/2 


1 - 


- 3 / 2 - 

'n? 


+ 0{e\ 

- 3 / 2-1 


0{e% 


(32) 

(33) 


where 'l'i(n) = neff(n,T)dT. Several comments are in 
order at this point. Firstly, the expressions of the two 
coefficients have exactly the expected form. When the 
parameter e goes to zero, Ci{n) goes to one and C 2 {n) 
goes to zero, i.e., we recover an unmodified spectrum. 
This is in complete agreement with the fact that, when 
e = rij(ri)/r]Q —> 0, the effective wavenumber UeS never 
penetrates the region where the WKB approximation is 
violated. This is in this sense that the parameter e is 


the relevant parameter to expand the spectrum in since 
it is clearly directly linked to the violation of the WKB 
approximation. Therefore, in the limit e —> 0, one re¬ 
covers the usual spectrum. Secondly, it is interesting 
to notice that the two terms proportional to e in Ci(n) 
and C 2 {n) are the same up to a global phase. To com¬ 
pute the phase in Eqs. (|2^) and (^, we have written 
fl(7/) = Iki)?!) -I- neff(n,T)dT. Since we only need to 
evaluate fl(±? 7 j), one can use the approximate equation 
















for rieff, see Eq. (|^), in the integral neff(n, r)dT. One 
finds 

n{±r]j) ='i>i±enr]oil -+C>(e^). (34) 

V ^■inf / 

Let us remark that the term in the ex¬ 

pression of C 2 (n) does not contribute to order e whereas 
the factor gives a contribution at this or¬ 

der. This contribution is exactly such that the two terms 
proportional to e in Eqs. (|2|), have the same ab¬ 
solute value. Thirdly, the factor 1 — is always 

positive since we restrict ourselves to the regime where 
Therefore, the above expressions of Ci{n) 
and (72 (n) are always well-defined. Finally, interestingly 


enough, the two coefficients Ci(n) and C' 2 (n) have a 
similar structure to one of the corresponding coefficients 
found in Ref. where another example of a modified 
spectrum has been exhibited. 

We are now in a position where we can calculate the 
spectrum explicitly. It is defined by the following equa¬ 
tion 


n^P(n) oc 



(35) 


We will evaluate the spectrum in the region where 
rj/rjo 1. In this case, airf) ~ €0 and UeS — n. Ex¬ 
panding everything in terms of the small parameter e, 
one obtains 


n^P{n) oc 


n'^ Vo 
211 '2^2 


n 





X sin 


2$(n) — 24'i(n) -I- 2n(?7 — 


where 77 + > is the time at which riefi — n and the 
phase $(n) is defined by $(n) = drrieff. Again the 
spectrum has the expected form. The leading order is the 
vacuum spectrum with spectral index 3 and the next-to- 
leading order (proportional to e) has the form of a com¬ 
plicated function of n times superimposed oscillations. 
We have plotted this complicated function (without the 
superimposed oscillations) in Fig. We can check that 



FIG. 5: The dependence of the matching time ? 7 j(n) on the 
wave number n for the values £0 = 5000, Ih = 10, = 

1, 1 5m I = 1, rn = 1. The solid curve is the exact result 
determined numerically, the dotted curve is obtained using 
the analytical approximations. 

the correction remains small and that, therefore, the ap¬ 





<f/ 



0 (e 2 ), 


(36) 


proximation employed in this article is consistent. On the 
other hand, when the wavenumber n approaches riinf, the 
approximation scheme breaks down but, precisely in this 
limit, the whole problem becomes pointless in agreement 
with the considerations above. We also see that the de¬ 
viation is minimum around n ~ 45 (for these values of 
the parameters) because, as already mentioned above, 
in this case the term a”/a almost never dominates. Al¬ 
though the present example is only a toy model, it is 
quite interesting to see that it bears some close resem¬ 
blance with the other example of a modified spectrum 
found in Ref. |Q. This suggests that the features previ¬ 
ously exhibited may be valid in general. 


We have thus reached our main goal: find an exam¬ 
ple where the initial conditions can be fixed naturally, 
where the frequency never becomes complex and where 
the final spectrum is modified. Note that in this exam¬ 
ple, the change in the spectrum comes about from an 
interplay between the modified dispersion relation factor 
rigjj and the factor a"/a which is responsible for Parker 
particle production |20 . For an unmodified linear dis¬ 


persion relation, the a"/a term is always negligible and 
there is no particle production. However, for our modi¬ 
fied dispersion relation, for a range of modes the Parker 
particle production term becomes important and leads 
to non-adiabatic evolution. The length of time during 
which the a"/a term dominates depends on the specifics 
of the dispersion relation, and hence the final spectrum 
will depend on these specifics. 












9 



FIG. 6: Amplitude of the next-to-leading order correction to 
the spectrum without the superimposed oscillations. Values 
chosen are £o = 5000, £b = 10, = 1, = 

1. When n approaches riinf the correction blows up and the 
approximation becomes meaningless. 

IV. DISCUSSION 

In this paper we have studied the dependence of the 
spectrum of a free scalar field in a bouncing Universe on 
trans-Planckian effects introduced via a modification of 
the free field dispersion relation. We have found that 
both the amplitude and the slope of the spectrum de¬ 
pend on the dispersion relation, assuming that the dis¬ 
persion relation leads to non-adiabatic evolution of the 
mode functions on trans-Planckian physical length scales. 
Such non-adiabatic evolution is possible without requir¬ 
ing the effective frequency to become imaginary. 

This result supports our earlier work which shows that 
the spectrum of cosmological fluctuations in inflation¬ 
ary cosmology can depend sensitively on trans-Planckian 
physics 01- Note that results supporting were 

recently also obtained by ^ (see also Ref. |^) in the 
context of a mode equation modified from the usual ones 
by taking into account effects coming from a short 
distance cutoff. For a recent paper on the effects of short- 
distance physics on the consistency relation for scalar and 
tensor fluctuations in inflationary cosmology see Ref. . 
The calculation we presented here is free of two of the 
possible objections against the earlier work. The first 


objection was that the use of the minimum energy den¬ 
sity state as initial state is not well defined on trans- 
Planckian scales if the dispersion relation differs dramat¬ 
ically from the linear one (which it has to in order to 
get non-adiabatic evolution). In our present work, the 
initial conditions are set in the low curvature region and 
on length scales larger than the cutoff length but smaller 
than the Hubble radius, where the choice of initial vac¬ 
uum state is well defined. The second objection con¬ 
cerned the use of dispersion relations for which the ef¬ 
fective frequency is imaginary in some time interval. No 
imaginary frequencies are used in the present model. 

The results obtained in this paper apply also both 
to bouncing Universe backgrounds described by the two 
other models mentioned in the Introduction. In the case 
of the first model (where the scale factor is given by an 
hyperbolic cosine), similar results would hold if we would 
match the local Minkowski vacuum (WKB vacuum) at 
some initial time and express the results in terms of the 
WKB vacuum state at the corresponding post-bounce 
time. However, in that case, the initial conditions are not 
easy to justify since they have to be set on length scales 
larger than the Hubble radius. Obviously, we could con¬ 
sider that model and let the scale factor make a further 
transition to an asymptotically radiation- dominated or 
asymptotically flat stage at very large initial and final 
times. In this case, it is very likely that results simi¬ 
lar to the ones we obtained here in model (^ would be 
obtained. 

Acknowledgments 

We would like to thank M. Lemoine for comments and 
careful reading of the manuscript. We acknowledge sup¬ 
port from the BROWN-CNRS University Accord which 
made possible the visit of J. M. to Brown during which 
most of the work on this project was done, and we are 
grateful to Herb Fried for his efforts to secure this Ac¬ 
cord. S. E. J. acknowledges financial support from CNPq. 
R. B. wishes to thank Bill Unruh for hospitality at the 
University of British Columbia during the time when this 
work was completed, and for many stimulating discus¬ 
sions. J. M. thanks the High Energy Group of Brown 
University for warm hospitality. The research was sup¬ 
ported in part by the U.S. Department of Energy under 
Contract DE-FG02-91ER40688, TASK A. 


[ 1 ] 

[ 2 ] 

[3] 

[4] 

[5] 


R. Branden berger and J. Marti n, Mod. Phys. Lett. A16, 
999 (2001), astro-ph/ 000543^ . 

J. Martin an d R. Brandenber ger. Phys. Rev. D63, 
123501 (2001), ^Ii^-th/000520^ . 

J. Niemeyer, Ph ys. Rev. D63, 123502 (2001), 

astro-ph/0005533|. 


W. Unruh, Phys. Rev. D51, 2827 (1995). 

S. Corley and T. Jacobson, Phys. Rev. D54, 1568 (1996); 


[ 6 ] 

[7] 

[ 8 ] 


S. Corley, Phys. Rev. D57, 6280 (1998). 

J. Kowalski-Glikrn an. Phys. Lett. B499, 1 (2001), 
astro-ph/0006250. 


J. Martin and R. Brandenberger, Proceedings of the 
Ninth Marcel Grossmann Meeting on General Relativity, 
edited by R. T. Jantzen, V. Gurz adyan and R. Ruffi ni. 
World Scientific, Singapore, 2002, astro-ph/0012031. 

J. Niemeyer and R. Parentani, Phys. Rev. D64, 101301 













10 


(2001), ^stro-ph/0101451 . 

[9] A. Sta robinskv, Pisma Z h. Eksp. Teor. Fiz. 73, 415 


(2001), istro-ph/010404; 


[10] M. Gas perini and G. Ven eziano. Astropart. Phys. 1, 317 
(1993), jiep-th /9211024 

[11] J. Khoury, B. A. Ovrut, P. J. St einhardt and N. Turok, 
Phys. Rev. 1)64. 123522 (200 1'). |hep-th/0103239, 

[12] D. H. Lyth, tiep-ph/0106153| . 


[13] R. Brandenberger and F. Finelli, iep-th/0109004 


[14] J. Khoury, B. A . Ovrut, P. J. Steinhardt and N. Turok, 
[hep-th/01 09050|. _ 


[15] J. Hwang, astro-ph/0109045 


[16] J. Martin. P. Pe ter. N. Pinto-Neto and D. J. Schwarz, 
|hep-th/0112128. 


[17] J. Audretsch and G. Schafer, Phys. Lett. 66A, 459 
(1978). 

[18] N. Birrell and P. Davies, ” Quantum fields in curved 
space!', Car nbridee Univ. Press . Cambridge, 1982. 

[19] T. Tanaka, astro-ph/0012431. 


[20] L. Parker, Phys. Rev. 183, 1057 (1969). 

[21] V. Mukhanov and R. Brandenberger, 
Phys. Rev. Lett. 68, 1969 (1992); R. Brandenberger, 
V. Mukhanov and A. Sornborger, Phys. Rev. D48, 1629 
(1993). 

[22] M. Brown and G. Dutton, Phys. Rev. D18, 4422 (1978). 

[23] L. Mersini, M. Bastero- Gil and P. Kanti, 


Phys. Rev. D64, 043508 (2001), hep-ph/0101210 


[24] M. Lemoine, M. Lubo, J. M ytin and J. P. Uzan, 


Phys. Rev. D65, 023510 (2002), hep-th/0109128 


[25] R. Easther, B. Greene, W. Kinney and G. Shin, 


Phys. Rev. D64, 103502 (2001), hep-th/0104102 


[26] A. Kert ipf and J. Niemey er. Phys. Rev. D64, 103501 
(2001), ^stro-ph/010325^. 

[27] A. Kempf, Phy s. Rev. D63, 083514 (2001), 


astro-ph/0009209 


L. Hui and W. H. Kinney, astro-ph/0109107 




































